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The attractor mechanism for supersymmetric (BPS) black holes was discov- 
ered in 1995 |T]: at the horizon of a supersymmetric black hole, the moduli are 
completely determined by the charges of the black hole, independent of their 
asymptotic values. In 2005, Sen showed that all extremal black holes, both su- 
persymmetric and non-supersymmetric (non-BPS), exhibit attractor behavior 
[2]: it is a result of the near-horizon geometry of extremal black holes, rather 
' than supersymmetry. Since then, non-BPS attractors have been a very active 

^ field of research (see for instance [3ll|51[ll[aiHl[S[Tni[IIl[Tl[Tl[Il[Tn]). 

In particular, a microstate counting for certain non-BPS black holes was pro- 
posed in [16] . Moreover, a new extension of topological string theory was 
J> suggested to generalize the Ooguri-Strominger-Vafa (OSV) formula so that it 

also applies to non-supersymmetric black holes [ITj. 

Both BPS and non-BPS attractor points are simply determined as the 
critical points of the black hole potential Vbh [HI [Z]- However, it is much 
easier to solve the full BPS attractor flow equations than to solve the non- 
' ^ BPS ones: the supersymmetry condition reduces the second-order equations 

of motion to first-order ones. Once the BPS attractor moduli are known in 
f"^ terms of D-brane charges, the full BPS attractor fiow can be generated via 

a harmonic function procedure, i.e., by replacing the charges in the attractor 
moduli with corresponding harmonic functions: 



> 

X 



^ tBPsix) = t*Bpsip' ^ H\x),qi ^ Hiix)) (1) 

In particular, when the harmonic functions [H^ (x), Hj{x)) are multi-centered, 
this procedure generates multi-centered BPS solutions [THj . 

The existence of multi-centered BPS bound states is crucial in under- 
standing the microscopic entropy counting of BPS black holes and the exact 
formulation of OSV formula [2QJ . One can imagine that a similarly important 
role could be played by multi-centered non-BPS solutions in understanding 
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non-BPS black holes microscopically. However, the multi-centered non-BPS 
attractor solutions have not been constructed until [21], on which this talk is 
based. In fact, even their existence has been in question. 

In the BPS case, the construction of multi-centered attractor solutions 
is a simple generalization of the full attractor flows of single-centered black 
holes: one needs simply to replace the single-centered harmonic functions in 
a single-centered BPS flow with multi-centered harmonic functions. However, 
the full attractor flow of a generic single-centered non-BPS black hole has not 
been solved analytically, due to the difficulty of solving second-order equations 
of motion. Ceresole et al. obtained an equivalent first-order equation for non- 
BPS attractors in terms of a "fake superpotential," but the fake superpotential 
can only be explicitly constructed for special charges and asymptotic moduli 
|22l I23j . Similarly, the harmonic function procedure was only shown to apply 
to a special subclass of non-BPS black holes, but has not been proven for 
generic cases [TT]. 

In this talk, we will develop a method of constructing generic black hole 
attractor solutions, both BPS and non-BPS, single-centered as well as multi- 
centered, in a large class of 4D Af — 2 supergravities coupled to vector- 
multiplets with cubic prepotentials. The method is applicable to models for 
which the 3D moduli spaces obtained via c*-map are symmetric coset spaces. 
All attractor solutions in such a 3D moduli space can be constructed alge- 
braically in a unified way. Then the 3D attractor solutions are mapped back 
into four dimensions to give 4D extremal black holes. 

The outline of the talk is as follows. Section 2 lays out the framework and 
presents our solution generating procedures; section 3 focuses on the theory 
of 4D Af = 2 supergravity coupled to one vector- multiplet, and shows in 
detail how to determine the attractor flow generators; section 4 then uses 
these generators to construct single-centered attractors, both BPS and non- 
BPS, and proves that generic non-BPS solutions cannot be generated via the 
harmonic function procedure; section 5 constructs multi-centered solutions, 
and shows the great contrasts between BPS and non-BPS ones. We end with 
a discussion on various future directions. 



2 Framework 

2.1 3D Moduli Space Msn 

The technique of studying stationary configurations of 4D supergravities by 
dimensionally reducing the 4D theories to 3D non-linear tr-models coupled to 
gravity was described in the pioneering work |24j . The 3D moduli space for 
4D AA = 2 supergravity coupled to ny vector-multiplets is well-studied, for 
example in [25l |26l |27l |28j . Here we briefly review the essential points. 
The bosonic part of the 4D action is: 
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S^-^ J d'^x^/-gi^) R - 2G,jdf A *4di^ - A G/ (2) 

where / = 0,l...ny, and Gj = {ReAf)i,jF'^ + {ImJ\f)ij * F-^ . For a theory 
endowed with a prepotential F{X), Afjj = Fjj + 2i ''^'^^ xlmF^x ^^ ' "^^^^re 
Fij — didjF{X) ■ We will consider generic stationary solutions, allowing 
non-zero angular momentum. The ansatz for the metric and gauge fields are: 

ds^ = -e2^(dt + u}f + e-^^g^fedx'^dx'' (3) 
A' ^ Al^{dt + uj) + (4) 

where g^,;, is the 3D space metric and bold fonts denote three-dimensional 
fields and operators. The variables are 3ny -I- 2 scalars {U,t^ ,t\ Aq}, and 
nv + 2 vectors {u:,A^}. 

The existence of a time- like isometry allows us to reduce the 4D theory to 
a 3D non-linear a-model on this isometry. Dualizing the vectors {u>,A^} to 
the scalars {a, Bj}, and renaming Aq as A^ , we arrive at the 3D Lagrangian, 
which is a non- linear cr-model minimally coupled to 3D gravity}^ 

^ = ^Vg(--R + a,0"9>"5™n) (5) 

Z K 

where (f>" are the 4(ny + 1) moduli fields {U,t^,t^,(j,A^ ,Bi}, and gmn is the 
metric of the 3D moduli space M3D1 whose line element is: 

ds^ = dU^ + ^e-'^"{da + A'dBi - BidA')^ + g,j{t, t)df ■ dP 



-e 



-^'^[{ImJ\f-y-'{dBj + AfjKdA'') ■ {dBj+AfjLdA'^)] (6) 



The resulting M3D is a para-quaternionic-Kahler manifold, with special 
holonomy 5*^(2, M) x 5*^(271^-1-2, M) [29J. It is the analytical continuation of the 
quaternionic-Kahler manifold with special holonomy USp{2,R) x USp{2nv + 
2, M.) studied in [26]. Thus the vielbein has two indices {a, A), transforming un- 
der Sp{2,M.) and 5*^(271^ -I- 2, M), respectively. The para-quaternionic vielbein 
is the analytical continuation of the quaternionic vielbein computed in .26j. 
This procedure is called the c*-map [53]) as it is the analytical continuation 
of the c-map in f25| [21] 

The isometrics of the Msd descends from the symmetry of the 4D system. 
In particular, the gauge symmetries in 4D give the shift isometries of Msd, 
whose associated conserved charges are: 

qjdr ^ Jj^i ^ P^i-BiPa, p^dr ^ Jbi = Pbi+A^ Pa, adr = = Pa 

(7) 



^ Note that the black hole potential term in 4D breaks down into kinetic terms of 
the 3D moduli, thus there is no potential term for the 3D moduli. 
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where the {Pa, Pa' , Pbj} are the momenta. Here r is the affine parameter 

defined as dr = — *^ smOdOdcf). {p^,qi) are the D-branc charges, and a the 
NUT charge. A non-zero a gives rise to closed time-hke curves, so we will set 
a = from now on. 

2.2 Attractor Flow Equations 

The E.O.M. of 3D gravity is Einstein's equation: 

Ra6 - ^g„6R = KTab = nidaCP"" dbC^^ Qmn ' Isabdc'P"' 9" ^ 9mn) (8) 

and the E.O.M. of the 3D moduli are the geodesic equations in Adao- 

V„V«<^" + Cp9a</.'"5«</)f = (9) 

It is not easy to solve a non-linear cr-niodel that couples to gravity. How- 
ever, the theory greatly simplifies when the 3D spatial slice is fiat: the dy- 
namics of the moduli are decoupled from that of 3D gravity: 

T„6 = = 5a0"96</)"w and dad'' ^ + r:},pda<P''' ^ = (10) 

In particular, a single-centered attractor flow then corresponds to a null 
geodesies in A^sij: ds'^ ~ dip^dcp^g^n = 0- 

The condition of the 3D spatial slice being fiat is guaranteed for BPS 
attractors, both single-centered and multi-centered, by supersymmetry. Fur- 
thermore, for single-centered attractors, both BPS and non-BPS, extremality 
condition ensures the fiatness of the 3D spatial slice. In this paper, we will 
impose this fiat 3D spatial slice condition on all multi-centered non-BPS at- 
tractors we are looking for. They correspond to the multi-centered solutions 
that are directly "assembled" by single-centered attractors, and have prop- 
erties similar to their single-centered constituents: they live in certain null 
totally geodesic sub-manifolds of A4sd- We will discuss the relaxation of this 
condition at the end of the paper. 

To summarize, the problem of finding 4D single-centered black hole attrac- 
tors can be translated into finding appropriate null geodesies in Mzd, and 
that of finding 4D multi-centered black hole bound states into finding cor- 
responding 3D multi-centered solutions living in certain null totally geodesic 
sub-manifold of M^d- 

The null geodesic that corresponds to a 4D black hole attractor is one that 
terminates at a point on the U —oo boundary and in the interior region 
with respect to all other coordinates of the moduli space A^3_d. However, it 
is difiicult to find such geodesies since a generic null geodesic fiows to the 
boundary of Msd- For BPS attractors, the termination of the null geodesic 
at its attractor point is guaranteed by the constraints imposed by supersym- 
metry. For non-BPS attractor, one need to find the constraints without the 
aid of supersymmetry. We will show that this can be done for models with 
A4:ii) that are symmetric coset spaces. Moreover, the method c;an be easily 
generalized to find the multi-centered non-BPS attractor solutions. 
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2.3 Models with A^sr) Being Symmetric Coset Spaces 

A homogeneous space is a manifold on which its isometry group G acts 
transitively. It is isomorphic to the coset space G/H, with G being the isom- 
etry group and H the isotropy group. For Msd ~ G/H, H is the maximal 
compact subgroup of G when one compactifies on a spatial isometry down to 
(1,2) space, or the analytical continuation of the maximal compact subgroup 
when one compactifies on the time isometry down to (0,3) space. 
The Lie algebra g has Cartan decomposition: g = h k where 

[h,h] = h [h,k]=k (11) 

When G is semi-simple, the coset space G/H is symmetric, meaning: 

[k,k] = h (12) 

The building block of the non-linear cr-model with symmetric coset space 
MsD as target space is the coset representative M, from which the left- 
invariant current is constructed: 

J ^ Ar^dM ^ Jk + Jh (13) 

where Jk is the projection of J onto the coset algebra k. The lagrangian 
density of the cr-model with target space G/H is then given by Jk as: 

/: = Tr(JkA*3Jk) (14) 

The symmetric coset space has the nice property that its geodesies M(t) 
are simply generated by exponentiation of the coset algebra k: 

M(t) = A-foe'^^/^ with fc e k (15) 

where Mq parameterizes the initial point of the geodesic, and the factor ^ in 
the exponent is for later convenience. A null geodesic corresponds to |fcp = 0. 
Therefore, in the symmetric coset space M3D, the problem of finding the 
null geodesies that terminate at attractor points is translated into finding the 
appropriate constraints on the null elements of the coset algebra k. 

The theories with 3D moduH spaces Man that are symmetric coset spaces 
include: £)-dimensional gravity toroidally compactified to four dimensions, all 
4D Af > 2 extended supergravities, and certain 4D M — 2 supergravities 
coupled to vector-multiplets with cubic prepotentials. The entropies in the 
last two classes are U-duality invariant. In this talk, we will focus on the last 
class. The discussion on the first class can be found in pi] . 

Parametrization of M-sd 

The symmetric coset space M^jj = G /H can be parameterized by exponenti- 
ation of the solvable subalgebra solv of g: 
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Al3D = G/H = e"°''' with g^h® solv (16) 

The solvable subalgebra solv is determined via Iwasawa decomposition of g. 
Being semi-simple, g has Iwasawa decomposition: g = h © a © n, where a 
is the maximal abelian subspace of k, and n the nilpotent subspace of the 
positive root space of a. The solvable subalgebra solv = a©n. Each point 
(/)'" in A^3D corresponds to a solvable element I^{(j>) = e*"'", thus the solvable 
elements can serve as coset representatives. 

We briefly explain how to extract the values of moduli from the coset 
representative M. Since M is defined up to the action of the isotropy group 
H, we need to construct from M an entity that encodes the values of moduli 
in an H-independent way. The symmetric matrix S defined as: 

S EE MSoM^ (17) 

has such a property, where Sq is the signature matrixj^ Moreover, as the 
isometry group G acts transitively on the space of matrices with signature 5*0, 
the space of possible 5* is the same as the symmetric coset space M^d = G/H. 
Therefore, we can read off the values of moduli from S in an H-independent 
way. 

The non- linear cr- model with target space Aiso can also be described in 
terms of S instead of M. First, the left-invariant current of S' is J5 = S~^dS, 
which is related to Jk by: 

Js = S-^dS = 2(S'oM^)-i Jk(5oAf^) (18) 

The lagrangian density in terms of S is thus C = ^Tr( J5 A*3 J5). The equation 
of motion is the conservation of current: 

V • J V • (S'-^VS') = (19) 

where we have dropped the subscript S* in Jg, since we will only be dealing 
with this current from now on. 

2.4 Example: Uy — 1 

In this talk, we will perform the explicit computation only for the simplest 
case: 4D J\f = 2 supergravity coupled to one vector-multiplet. The generaliza- 
tion to generic ny is straightforward. The 3D moduli space M^d for ny — 1 
is an eight-dimensional quaternionic kahler manifold, with special holonomy 
5*^(2, M) X 5*^(4,11). Computing the killing symmetries of the metric (|6| with 
nv = 1 shows that it is a coset space G2{2)/{SL{2,R) x S'L(2,M))[3] Figure[l] 
shows the root diagram of 6*2(2) in its Cartan decomposition. The six roots on 

^ In all systems considered in the present work, the isotropy group H is the maximal 
orthogonal subgroup of G: HS^H^ = So, for any G H. Therefore, S is invariant 
under the H-action M ^ MH. 

^ Other work on this coset space has appeared recently, including |31l 1321 133] . 
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Fig. 1. Root Diagram of G2(2) in Cartan Decomposition. 
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the horizontal and vertical axes {L^, L^, ij, L^} generate the isotropy sub- 
group H = SL{2, R)fi X SL{2, M.)^. The two vertical columns of eight roots QaA 
generate the coset algebra k, with index a labeling a spin-1/2 representation 
of SL{2,M.)h and index A a spin-3/2 representation of SL{2,M.)y. 

The Iwasawa decomposition, g = h solv with solv — a n, is 
shown in Figure [2] The two Cartan generators {u, y} form a, while n is 
spanned by {x, cr, A", A^, Bi, Bo}, {u, y} generates the rescahng of {u,y} 
where u = e 

{x,a,A°,A\B,,Bo}m 

The moduli space Msd can be parameterized by solvable elements: 



^'^ and {x, cr, A", A^, Bi, Bq} generates the translation of 



J\nu)u/2+(lny)y xx+A' a' + B,Bi+arT 



(20) 



The symmetric matrix S can then be expressed in terms of the eight moduli 



(21) 



which shows how to extract the values of moduli from S even when S is not 
constructed from the solvable elements, since it is invariant under H-action. 
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Fig. 2 

{u, y, X, (T, A*^, A^, Bi, Bo} generates the solvable subgroup 5*0/11. 



Root Diagram of G 



2(2) 



A" 



in Iwasawa Decomposition. 



3 Generators of Attractor Flows 



In this section, we will solve 3D attractor flow generators k as in (15 1. We 
will prove that extremality condition ensures that they are nilpotent elements 
of the coset algebra k. In particular, for ny = 1, both BPS and non-BPS 
generators are third-degree nilpotent. However, despite this common feature, 
ksps and fcjvs differ in many aspects. 



3.1 Construction of Attractor Flow Generators 



Construction of ksps 

Since the 4D BPS attractor solutions are already known, one can easily obtain 
the BPS flow generator ksps in the 3D moduli space A^sd. 

The generator ksps can be expanded by coset elements aaA as ksps — 
QaAC"^, where C"^ are conserved along the flow. On the other hand, since 
the conserved currents in the homogeneous space are constructed by project- 
ing the one-form valued Lie algebra -dg onto k, a procedure that also gives 
the vielbein: Jk = g~^dg\]i = UaAV"^, the vielbein V"^ are also conserved 

along the flow: (^'^0") = 0. Since both the expansion coefficients C"^ 
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and the viclbein V"^ transform as (2,4) of 5L(2,M)/, x (5L(2,M)„ and are 
conserved along the flow, they are related by: 

(jo^A ^ v^^<i)" (22) 

up to an overall scaling factor. 

In terms of the vielbein V"^, the supersymmetry condition that gives the 



BPS attractors is: V""^ = z"T/^ [2911301133]. Using (|22|, we conclude that the 
3D BPS flow generator ksps has the expansion 

ksps = aa,Az"C^ (23) 

A 4D supersymmetric black hole is labeled by four D-brane charges 
{p^,p^, gi, go)- A 3D attractor flow generator ksps has flve parameters {C^, z}. 
As will be shown later, z drops off in the final solutions of BPS attractor flows, 
under the zero NUT charge condition. Thus the geodesies generated by ksps 
are indeed in a four-parameter family. 

ksps can be obtained by a twisting procedure as follows. First, define 
a kgpg which is spanned by the four coset generators with positive charges 
under S'L(2,M),,: 

k%ps = aiAC^ (24) 
then, conjugate kgpg with lowering operator LJ^: 

kBPS ^ e-'''^^ k^pse'"^^ (25) 

Using properties of kppg, it is easy to check that kgpg is null: 

|fcBPs|'=0 (26) 

More importantly, ksps is found to be third-degree nilpotent: 

k%ps = (27) 

A natural question then arises: Is the nilpotency condition of ksps a result of 
supersymmetry or extremality? If latter, we can use the nilpotency condition 
as a constraint to solve for the non-BPS attractor generators k^p. We will 
prove that this is indeed the case. 



Extremality implies nilpotency of flow generators 

We will now prove that all attractor fiow generators, both BPS and Non-BPS, 
are nilpotent elements in the coset algebra k. It is a result of the near-horizon 
geometry of extremal black holes. 

The near- horizon geometry of a iD attractor is AdS2 x S'^, i.e. 

as r ^ oo (28) 
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As the flow goes to the near- horizon, i.e., as u ~ e^^ 0, the solvable element 
M — e(^"")"/2+- - is a polynomial function of r: 

M{t) - u-^/^ - (29) 

where — £ is the lowest eigenvalue of u. 

On the other hand, since the geodesic flow is generated by € k via 
M(r) = Moe'=^/^ M(r) is an exponential function of r. To reconcile the two 
statements, the attractor flow generator k must be nilpotent: 

fc^+i = (30) 

where the value of i depends on the particular moduli space under consid- 
eration. In G'2(2)/<S'L(2, K)^, by looking at the weights of the fundamental 
representation, we see that t = 2, thus 

k^ = (31) 

The nilpotency condition of the flow generators also automatically guar- 
antees that they are null: 

P = =^ [k^f =0 => Tr(fc2) = (32) 



Construction of fcATB 

To construct non-BPS attractor flows, one needs to flnd third-degree nilpotent 
elements in the coset algebra k that are distinct from the BPS ones. In the 
real 6*2(2) /'S'i(2,M)^, there are two third-degree nilpotent orbits in total ,35] . 

We have shown that fcsps ~ e^^''^ k^p^e^^i^ , with kgpg spanned by the 
four generators with positive charge under SL{2,R)h- 

Since there are only two S'L(2,]R)'s inside H, a natural guess for k^s is 
that it can be constructed by the same twisting procedure with SL{2,R)h 
replaced by S'i(2,M)„: 

kMB^e-'^^^k^ge'^^ with k%B=aa,aC'"', a,a = l,2 (33) 

where fc^^ is spanned by the four generators with positive charge under 
SL(2,M)„. 

Using properties of fcjvs' easily show that kj^s defined above is 

indeed third-degree nilpotent: 

k%B = (34) 



That is, k^B defined in (33) generates non-BPS attractor flows in Adso- 

A 4D non-BPS extremal black hole is labeled by four D-brane charges 
{p^,p^,qi,qo)- Similar to the BPS case, the 3D attractor flow generator fcjvs 
has five parameters {C°"^,z}. As will be shown later, z can be determined 
in terms of {C""} using the zero NUT charge condition, thus the geodesies 
generated by k]\rB are also in a four-parameter family. 



Non-Supersymmetric Attractors in Symmetric Coset Spaces 



11 



3.2 Properties of Attractor Flow Generators 

We choose the representation of 6*2(2) group to be the symmetric 7x7 
matrices that preserve a non-degenerate three- form Wijk such that rjis = 
WijkW stuWmnof--'''^^™"'" IS a metric with signature (4, 3) and normaHzed as 
rf' = 1. We decompose 7 as 3 3 1 of S'L(3,M) and choose the non-zero 
components ofw, 3A3A3, 3A3A3 and 3 (g) 3 (g) 1, as 

w = dxi A dx2 A -t- dy^ A dy'^ A dy^ \=dxa A dy"" A dz (35) 

v2 

which gives 77 = dxady'^ — dz^. Written exphcitly, an element of 6*2(2) Lie 
algebra is 

/ All ^nnkv'' V2w^, \ 
g = e'2jifey;^ _A^2^ _^y^2 (36) 

V-%/2wJi V2wj^ / 

Here A is a traceless 3x3 matrix. The signature matrix Sq is thus Diag [1,— 1,— 1,1,— 1,— 1,1]. 

The real 6*2(2) group has two third-degree nilpotent orbits. In both orbits, 
fc^ is of rank two and has Jordan form with two blocks of size 3. Thus can 
be written as 

/C^ = ^ Vav'[cabSo (37) 
a,b=l,2 

with Va null and orthogonal to each other: Va ■ Vb = 'v^S^Vb — 0, and Cab 
depends on the particular choice of k. Therefore, k can be expressed as: 

fc = ^ {Vawl + Wavl)So (38) 

a=l,2 

where each Wa is orthogonal to both Va'- Wa ■ Vb = 0, and Wa satisfies Wa • 
Wb — Cab- Next we solve for Va and Wa for ksps arid k^B and compare their 
properties. 



Properties of ksps 

The null space of k^ is five-dimensional, with Vg s panning its two-dimensional 
complement. For ksps, ^a^^ arid w^^^ in (38l are solved in terms of 
and z. 

In basis (36 1, from inspection of k^pg, we find that v^^^ can always be 
chosen to have the form0 

^^r^'^ (^1,-^1^1,0) vr"" = {-V2,mV2,V2) (39) 



There are some freedom on the choice of (va,Wa): a rotational freedom: 
iva^Wa) —* (RabVb, RabWb) with R orthogonal; and a rescaling freedom: 

{Va,'Wa) -+ {Var,Wa/r). 



12 Wei Li 



where rji is a 3D signature matrix 77 = Diag[l, —1,-1], and Va are two three- 
vectors satisfying 

Vi-Vi^O Vi-V2^0 V2-V2 = -I (40) 

We drop the superscript "BPS" for Va here since, as will be shown later, v^^ 
can also be written in terms of Va, though in a slightly different form. Note 
that for ksps, V2 is defined up to a shift of Vi: V2 ^ V2 — CV1, since any linear 
combination of Va^^ forms a new set of vj^^^ . 

Written in twistor representation are given by the twistors z and u as 

VI''' = 2z''z^ V:^^ = + z^u'' (41) 

where we have used the rescaling freedom to set z^v? — z^v}- — 1. Note that 
for kspSi the twistor u is arbitrary, due to the shift freedom of V2. 
The condition w^^^ ■ v^^^ = dictates that w^^^ has the form: 

= (T4^f ViWr', 0) w^P^ = (M^2^^^, r/iM^f 0) (42) 

with Wa^^ solved as: 

{W^^^, Wi^^)""^ = (P'^^-'u^^P'^^^z^) (43) 

where the totally symmetric P"P'> is defined in terms of as 

pl" = P"2 ^ ^2 pl22 ^ (j3 p222 ^ ^4 ^^^-^ 

In summary, v^^^ span a one-dimensional space (since u is arbitrary) and 
w^^^ span a four-dimensional space. 

3.3 Properties of k^s 

{v^^ ,Wa^) are solved in terms of {C"", z}. The forms of v^^ are only slightly 
different from those of v^^^: 



,NB 



{Vi,mVi,o) <^ - {V2, -mV2, V2) (45) 



where Va are the same three- vectors given in (41 1, with one major difference: 



the twistor u is no longer arbitrary, but is determined by C"" as 

^ With the inner product of two three- vectors defined as Va ■ Vb = V^rjiVb, the 
twistor representation of a three-vector V = {x, y, z) can be chosen as 



X + y z 
z X — y 
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since the V2 in v^^ no longer has the shift freedom. 

The forms of w^^ are also only slightly different from the BPS ones ( 42 1 : 

^NB ^ ^wr,-mwr,o) <^ = {wi'^,mwr,o) (47) 

with W^^ solved in terms of {C"°, z,u} as: 



Since the value of u imposes an extra constraint on the vectors w^^ via ( 46 1 , 
w^^ span a three-dimensional space instead of a four-dimensional one as in 
the BPS case (43 1. In summary, in contrast to the BPS case, v^^ span a 
two-dimensional space and w^^ span a three-dimensional one. 



4 Single-centered Attractor Flows 

Having solved the attractor flow generators for both BPS and non-BPS case, 
we are ready to construct single-centered attractor flows. A geodesic starting 
from arbitrary asymptotic moduli is given by M{t) = Afoe'''^/2, which gives 
the flow of S as S{t) = Mqc'''^ SqMq , which in turn can be written as S{t) = 
e^^'^^So, where K{t) is a matrix function. From now on, we use capital K 
to denote the matrix function which we exponentiate to generate attractor 
solutions. 

The current of S is 

J = S-^\7S ^ So (s/K+[S/K,K] + ^[[S/K,K],K] + ---^ So (50) 

The equation of motion is the conservation of currents: V • (S'~^V5') — 0, 
which is solved by K{t) being harmonic: 

V^X(t) = => K{T) = kT + g (51) 

g parameterizes the asymptotic moduli. Using the H-action, we can adjust 
g such that g € k, and g has the same properties as the flow generator k, 
namely, g^ = and g^ is of rank two. Therefore, for single-centered flow 
given by S{t) = e^^'^^So, the harmonic matrix function K{t) has the same 
properties as the flow generator k: 

K^{t)^0 and j.^^]^ (52) 



To find the harmonic K{t) that satisfies the constraints (52), recall that 
the constraints dictate K{t) to have the form: 

K{t) - MT)Wa{Tf + Wa{T)Va{Tf)So (53) 



14 Wei Li 



with Va{T) being null and Wa{T) orthogonal to Wfc(r) for all r. Then the con- 



straints (52 1 can simply be solved by choosing t'a(r) to be the constant null 
vectors Va{T) = Va and Wa^r) to be harmonic vectors which are everywhere 
orthogonal to w^: 

Wair) = WaT + rUa with Wa ■ Vb ^ ma ■ Vb ^ (54) 

The two 7-vectors WaS contain the information of the black hole charges, and 
the two 7-vectors Wa's contain that of asymptotic moduli. 

To summarize, the single-centered attractor flow starting from an arbitrary 
asymptotic moduli is generated by S{t) = e^^'^-'S'o, with harmonic matrix 
function K{t) = kr + g where 

k= [vaWa + WaVa]So and g= ^ [vainl + mavJ]SQ (55) 

a=l,2 a=l,2 

Since k and g share the same set of null vectors Va and both Wa and TOq are 
orthogonal to Vb, g has the same form as that of flow generator fc, namely: 

gsps = a^Az"G^ gNB - e"-^^' (a„„G"'^)e^^^' (56) 

which guarantees that g is also third-degree nilpotent. Moreover, that g and 
k have the same form implies = 0, thus the current is reduced to 

from which we can solve Va and Wa in terms of charges and asymptotic moduli. 

Now that we are able to construct arbitrary attractor flows in the 3D 
moduli space, we can lift them to the 4D black hole attractor solution. First, 



in representation given by ( 36 ) , the 4D moduli t = x + iy can be extracted 



from the symmetric matrix S via: 



(59) 



^ 'S'35(t) S33{t)S55{t) - 535(r)^ 

^^^^ = -s^) ^''^ 

and u — e^^ via: 

1 

7^33(r)555(r)-5^^ 

Since both k and g are third-degree nilpotent, S{t) is a quadratic function of 
r. Moreover, since g has the same form as fc, S{t) is composed of harmonic 
functions of r: H^{t) = C^t + G^ for BPS attractors and i7"°(r) = C°"'t + 
Qaa non-BPS attractorsj^] Generic single-centered attractor flows with 
arbitrary charges and asymptotic moduli can thus be generated. The attractor 



^ Space prohibits listing the rather lengthy result of S(t), readers can consult eq. 
(6.2) and eq. (6.3) of [2T1 for BPS attractors, and eq. (6.28) for non-BPS ones. 
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moduli are read off from S{t) with t ~f oo, and asymptotic moduli with 
0. 

The D-brane charges can be read off from the charge matrix defined as 
Q = J V • J. The AD gauge currents sit in the current J = S^-^\7S as: 

- 2 - 

(Jai, J51, J72, Ji2, J32) = {V2Jbo,-V2Jbi, - Jai, V2Jao, -2Ja) (60) 

Therefore Q relates to the D-brane charge {p^ , , qi , qo) and the vanishing 
NUT charge a by 

(Q3i,Q5i,Q72,Qi2) = (y2/,-y2p\^<7i,V2go) Q32 = -2a = (61) 



4.1 Single- Centered BPS Attractor Flows 

As an example, a single-centered BPS black hole constructed by lifting the 
attractor solution in Ai^o is shown in Figure |3] It has D-brane charges 
{p'^,p^,qi,qQ) — (5,2,7,-3). The four flows, starting from different asymp- 
totic moduli, terminate at the attractor point {xgpg^ygpg) with different 
tangent directions. The reason is that the mass matrix of the black hole po- 
tential Vbh at the BPS critical point has two identical eigenvalues, thus there 
is no preferred direction for the geodesies to flow to the attractor point. 

We now discuss in detail how to determine ksps and gsps for given 
charges and asymptotic moduli. There are nine parameters in ksps and gsps'- 
{C^ ^G^ , z}, since the twistor u is arbitrary. On the other hand, there are 
eight constraints in a given attractor flow: four D-brane charges {p^,qi), the 
vanishing NUT charge a, and the asymptotic moduli (xg, j/oi ''^o)!^ We will 
use these eight constraints to fix and in ksps and gspsi leaving the 
twistor z unfixed. 



Integrating the current (57 1 for BPS case produces five coupled equations: 



Qsps = Soiksps + 2i^BPS,gBPs])Sa (62) 

where [kBPS,gBPs] = (C, G)6>, with (C,G) = C^G^-3C'^G^+3C^G'^-C^G\ 
and e = -j^e-^'^^L+e^^'f^ . 

In order to show that the BPS fiow can be expressed in terms of harmonic 
functions: H{t) — Qt + h, with Q = {p^ , qj) and h = {h^ , hj), we will solve 
gsps in terms of h instead of (ccq, yoj mq). h relates to the asymptotic moduli 

by 

{xQ,yo, uq)bps = {x, y, u)'bps{Q h) (63) 
and there is one extra degree of freedom to be fixed later. 



The cisymptotic value of u can be fixed to an arbitrary value by a rescaling of 
time and the radial distance. We will set uo = 1- 
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First, for later convenience, we separate from gsps a piece that has the 
same dependence on {h, z) as ksps on (Q,z): 

9bps 9BPSM + A with gBPs.h = kspsiQ ^ h, z) (64) 

that is, gsps.h — o-aAZ^Gh ^i^h = C^{Q — > h). We can use the unfixed 
degree of freedom in h to set (C, Gh) — 0, so that (|62]) simpHfies into 



Qbps = SoikBPs + ^[kBPS,A])So (65) 



A can then be determined using the three constraints from ( 63 1 and the zero 
NUT charge condition in (lesj): A = a^Az'^E^ with = -E^^ 
E'^ = —E'^ = 1^5" • '^^^ form of A will ensure that the twistor z drops off in 
the final attractor flow solution written in terms of Q and h. 
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The remaining four conditions in the coupled equations ([65| determine 
as functions of D-brane charges and the twistor z: — C^C}, z) j^Then 
are given by = G^{Q h,z). The product {C^,G^) is proportional to 
the symplectic product of {p^,qi) and {h^,hj): 

2 

(G^, Gu) = ^ <Q,h> where <Q,h>= p^ho+p^hi-qih^-qoh° 

I + 

(66) 

The condition {G^,G^) = is then the integrability condition on h: < 
Q,h>^0. 

BPS attractor flows in terms of {p^ , qi) and {h^ , hj) are obtained by substi- 
tuting solutions of G^{Q, z) and G^{h, z) into the flow of S{t). The attractor 
moduli are determined by the charges as: 

, Ao + P^f , _ v/J4(p°,p\f,go) 

^^^^^ 2[(pi)2+p0^] yBPS- 2[(pi)2+p0^] (^7) 
where J4(p",p^, gi, go) is the quartic i?7(7) invariant: 
J4(p°,p\(?i,(7o) = 3(pigi)2-6(p°(7o)(p'gi)-(p°'Zo)'-4(pi)^go+V(<Zi)^ (68) 

thus J4(p°,p^, go) is the discriminant of charges. Charges with positive 
(negative) Ji{p'^,p^, ^j^o) form a BPS (non-BPS) black hole. The attractor 
value of = l/-\/^4(p°,P^, go)- The constraint on h from mo = 1 is 



then Ji{h^ , h^, /lo) = 1. The attractor moduli (67 1 match those from Type 



II string compactified on diagonal T^, with qi 

Now we will prove that the BPS attractor flows constructed above can 
indeed be generated by the "naive" harmonic function procedure, namely, by 
replacing charges Q in the attractor moduli with the corresponding harmonic 
functions Qt + h. First, using the properties of A, the flow of t — x + iy can 
be generated from the attractor moduli by replacing ps with the harmonic 
function kspsT + gspsX- 

tBPs{T) = t*BPs{^BPS kr + gBPS,h) (69) 

Then, since ksps and gBPS,h share the same twistor z, this is equivalent to 
replacing G^ with harmonic functions H^{t) = C'^r + G^ while keeping the 
twistor z fixed: 

tBPsir)^t*sps{C''-^C^T + Gtz) (70) 

Finally, since G^ is linear in Q and G^ linear in h, and since z drops off after 
plugging in the solutions G^{Q, z) and G'^{h, z), we conclude that the flow of 
tBPsiT) is given by replacing the charges Q in the attractor moduli with the 
corresponding harmonic functions Qt + h: 

tBPs{T)=t*Bps{Q^QT + h) (71) 
* See eq. (6.18) of [2T] for the full solutions. 
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4.2 Single- Centered Non-BPS Attractor Flows 

A non-BPS attractor flow with generic charges and asymptotic moduli can 
be generated using the method detailed earlier. Figure |4] shows an example 
of non-BPS attractor flow with charges {p'^,p^, qi, qq) = (5, 2, 7, 3). Note that 
J4(5, 2, 7/3, 3) < 0, so this is indeed a non-BPS black hole. 



y 




Fig. 4. Non-BPS flow with charges {p° , , qi , qo) = (5,2,7,3) and attractor point 
{x*,y*) = (—0.323385,0.580375). The initial points of each flow are given by: {xi = 
0.539624,7/1 = 5.461135), (a;2 = 1.67984,1/2 = 0.518725), (a;3 = -0.432811,7/3 = 
0.289493), (x4 = 1.28447, y4 = 1.49815), (xs = -0.499491,1/5 = 0.181744). 



Unlike the BPS attractor flows, all non-BPS flows starting from different 
asymptotic moduli reach the attractor point with the same tangent direction. 
The reason is, unlike the BPS case, the mass matrix of the black-hole potential 
Vbh at a non-BPS critical point has two different eigenvalues. The common 
tangent direction for the non-BPS flows corresponds to the eigenvector asso- 
ciated with the smaller mass. 
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Now we discuss how to determine b and g^B for given D-brane charges 
and asymptotic moduli. Unlike the BPS case, there are only eight parameters 
in kjsiB and gNB- the two twistors {z, u} and {C"°, G""} under the constraints 
u — — ^12- On the other hand, there are still eight constraints in a given 
non-BPS attractor flow as in the BPS case. Therefore, while ksps and gBPS 
can parameterize black holes with arbitrary (p^ , qj) and {xq, yo) while leaving 
{z, u} free, all the parameters in kj^g and g^Bi including {z, u}, will be fixed. 

Another major difference from the BPS case is that 

[kNB,gNB]^0 (72) 



guaranteed by the form of v^^ and w^^. Thus the charge equation (62) 
becomes simply 

Qnb = So{kNB)So (73) 

Unlike the BPS case, g^B does not enter the charge equations, thus cannot 
be used to eliminate the dependence on the twistor z. The three degrees 
of freedom in gpjB are simply fixed by the asymptotic moduli (xo^uq) and 
Up = 1, without invoking the zero NUT charge condition. The four D-brane 
charges equations in (73) determine — C"''((5, z)j^ which then fixes 



u 



via u — ^^712- Finally, the zero NUT charge condition imposes a degree-six 
equation on twistor z: 



P 



■°z^+6p^z^^{'3p"+Aqi)z^-4{3p^-2qa)z^ + {3p"+Aqi)z^+6p^z-p° ^0 (74) 



Similar to the BPS case, the full non-BPS attractor flow can be gener- 
ated from the attractor moduli by replacing with the harmonic function 



H°'°'{t) = C^^T + G"", while keeping z fixed as in (70 1 



^^^B(T) = ^^s(C"'^^C"V + G"^z) (75) 

However, there are two important differences. First, the harmonic functions 
H°"^ have to satisfy the constraint p°] 

Second, unlike the BPS flow, a generic non-BPS flow cannot be given by 
the "naive" harmonic function procedure: 

tNB (r) ^ tl^B {Q^Qt + h) (77) 



See eq. (6.35) of [21 for the full solution. 

This does not impose any constraint on the allowed asymptotic moduli since there 
are still three degrees of freedom in C^" to account for (a;o, j/o, '"o)- We will see 
later that its multi-centered counterpart helps impose a stringent constraint on 
the allowed D-brane charges in multi-centered non-BPS solutions. 
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The reason is that the twistor z in a non-BPS solution is no longer free as in 



the BPS case, but is determined in terms of D-brane charges via (74). Thus 
replacing Q with Qr+h, for generic Q and h, would not leave z invariant. That 
is, replacing C°"^ in the attractor moduli with harmonic functions H°'°'{t) is 
not equivalent to replacing the charges Q with H ~ Qt + as in the BPS 



case ( 71 ) 



It is interesting to find the subset of non-BPS single-centered fiows that 
can be constructed via the "naive" harmonic function procedure. The ny — 1 
system can be considered as the STU model with the three moduli {S,T,U) 
identified. Since the STU model has an S'L(2,Z)^ duality symmetry at the 
level of E.O.M., the riy = 1 system has an SL{2, Z) duality symmetry coming 

from identifying these three 5*27(2, Z)'s, namely, = " " '^^ 



c d J \cdj \cd 
with ad — be ^ 1 . The modulus t ^ x + iy transforms as t ^ Ft ~ ^^|qt^ , and 
the transformation on the charges is given by |34j . 

Given an arbitrary charge Q, there exists a transformation Fq such that 
Q = FqQ/^q for some D4-D0 charge system Q40 — (0,p^,0, go)- The solution 



of (74 1 with charge Q — FQQii^ has a root z ~ a±^/a-+c^ _^ independent of Q40. 
Thus for arbitrary hio = (0, /i^, 0, /iq), replacing Q with Qt -\- FqU/^q would 
leave the twistor z invariant. We thus conclude that the non-BPS single- 
centered attractor fiows that can be generated from their attractor moduli 
via the "naive" harmonic function procedure are only those with (Q, h) being 
the image of a single transformation F on a, D4-D0 system ((54o, ^■40): 

tNB (r) = t*j,B [rQio rQioT + rh^o) (78) 



5 Multi-Centered Attractor Flows 

Similar to the single-centered attractor solutions, the multi-centered ones are 
constructed by exponentiating harmonic matrix functions K{x): 

S{x) = e^^^'^o (79) 

Recall that for single-centered attractors, using the H-action on g, K{t) — 
kr + g can be adjusted to have the same properties as the fiow generator k 



as in (52 1. For BPS multi-centered solutions, supersymmetry guarantees that 



the matrix function K(x) also has the same properties as the generator k: 

K^{x) = and K'^{x) rank two (80) 

We will impose these constraints on all non-BPS multi-centered solutions as 
well, since presently we are more interested in the multi-centered solutions that 
are "assembled" by individual single-centered attractors and thus have similar 
properties to their single-centered constituents. It is certainly interesting to 
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see if there exist non-BPS multi-centered solutions with K{x) not sharing the 
constraints (80 1 satisfied by the flow generator fcjvs- 



The harmonic matrix function K(x) satisfying all the above constraints is 
solved to be: 

where 

'^[va{wa)J' + {wa)^v'^]So and ^ [vaml+mav'^]So (82) 

a=l,2 a=l,2 

with Va being the same two constant null vectors in single-centered fc, and the 
7- vectors {wa)i contain the information of the D-brane charges of center-i, 
and the two 7- vectors ma's contain that of asymptotic moduli. Both {wa)i 
and iria are orthogonal to Vf,. Since Va only depends on the twistor {z,u}, 
and Wa are linear in or C"", the above generating procedure is equivalent 
to replacing and C"" with the multi-centered harmonic functions H^{x) 
and H"°'{x) while keeping the twistor {z,u} fixed. 

5.1 Multi- Centered BPS Attractors 

Using Qi to denote the charge matrix of center-?, we have 5N coupled equa- 
tions from = 4^ /i V • J: 



Qbps^i = ^o{kBPS,i + 7)V^bps,i,9bps\ + 9 2^ f-^: — 1 (83) 

j 

We now show in detail how to determine kspsA and qbps for given charges 
and asymptotic moduli using equation (83). There are 4(Af+l) + l parameters 
in ksps.i find qbps'- {C'i^j G"*, z}, since the twistor u is arbitrary. Different 
from the single-centered BPS case, there are also ?>N — 3 degrees of freedom 



from the positions of centers on L.H.S. of (83 1. On the other hand, there 
are 5A^ -I- 3 constraints in a given BPS multi-centered attractor: AN D-brane 
charges {pi,qi.i), N vanishing NUT charges, the asymptotic moduli {xo,yo) 
and uq = 1. We will use these 5A^-|-3 constraints to fix the 4(A^-|-1) parameters 
{Cf, G^} in kBPS,i and gBPS: and impose iV— 1 constraints on the distances 
between the N centers, while leaving the twistor z free. 

First, integrating V • J over the sphere at the infinity gives the sum of 
the above N matrix equations: Q^ps ~ ^oi^^BPS + ^[^bps^ 9bps])So, which 
is the same as the charge equation for a single-center attractor with charge 
Qbps- '^^^^ determines g to he g — + A, same as the single-centered case 



as in (64 1, using the three asymptotic moduli (a;o,yoi'"o) and the constraint 
of zero total NUT charge. The h's are fixed by the asymptotic moduli and the 
integrability condition < Q*^pg,h >— 0. 



22 Wei Li 



It is easy to see that the solutions of Cf are simply given by the single- 
centered solutions = C^{Q,z) with Q replaced by Qj. Thus the flow 
generator of each center kBPs,i (given by kBPss = o-aAZ^Cf) satisfies 

Qbpsa = So{kBPs,i+ ^[kBPS,i,^])So (84) 
which is the multi-centered generalization of the single-centered condition 



Using the solutions of kBPS.i and gBPSi the charge equations (83) become 

QbPS.^S, I kBPS.:+\[kBPS,,A] + {< QBPS,.h> + ^ ^ QbPS,,QbPS,j > 
Z ^ \Xi Xj\ 

(85) 

from which we subtract ( 84 1 to produce the integrability condition 

^ ^ ^ I QBPS,i,QBPS,j > „ fQa\ 

< QBPS.i,h > + ■ ■ — - — = (86) 



The sum of the N equations in the integrability condition ( 86 1 reproduces 
the constraint on h: < Q*^pg,h >= 0. Thus the remaining N — 1 equations 
impose A^ — 1 constraints on the relative positions between the N centers. The 
angular momentum J, defined via lui — ie^j^J^^ as r — > oo, is non-zero: 

-J =\T.^^;^\^Qbps,,Qbps,,) (87) 

Thus we have shown that our multi-centered BPS attractor solutions repro- 
duce those found in [19j. Same arguments as in the single-centered BPS case 
shows that multi-centered BPS attractors can be generated by replacing the 
charges in the attractor moduli with corresponding multi-centered harmonic 
functions: 

tBPs{x) = t*Bps{QBPS X! I ^'^'j + ^) 



5.2 Multi- Centered Non-BPS Attractors 

A multi-centered non-BPS attractor has 3(iV-|-l)+2 parameters inside its non- 



BPS generators {/catb,*, Stvs}: {Cf^jG""} under the constraint (89l plus two 
twistors {z, u}. Given {fcjvs i, qnb} in terms of {Cf° , G"°, u], the non-BPS 
multi-centered solution is the same as the single-centered one with H°"^{t) 
replaced by multi-centered harmonic functions H"°'{x) = | -I- 

satisfying the constraint 

HPix) CP G12 ^^^^ 
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However, the process of determining k^B,! and Qnb in terms of charges and 
asymptotic moduli for a non-BPS multi-centered attractor is very different 
from its BPS counterpart. 

The reason is that the charge equations for a non-BPS multi-centered 
solution simplifies a great deal since 

[fcjvs,*, fcATBj] = and [fcwB.i, g^s] = (90) 

guaranteed by the forms of {w^^)i and m^^. Therefore, the 5N equations 



(83) decouple into N sets of 5 coupled equations: 

Qws,i = '5'o(fcjvB,i)S'o (91) 
As in the single-centered non-BPS case, g^B does not enter the charge equa- 



tions (91), and its three degrees of freedom can be completely fixed by the 
given asymptotic moduli (a;o,yo) and mq = 1 without using the zero NUT 
charge condition. More importantly, unlike BPS multi-centered solutions, the 



positions of centers Xi do not appear in the charge equations (91 ), thus receive 
no constraint: all centers are free. Finally, since we are using the remaining 
3A^ -I- 2 parameters {C"", z, u} to parameterize a A^-centered attractor solu- 



tion under 5A'^ constraints coming from charge equations (91), there need to 
be 2N — 2 constraints imposed on the D-brane charges. 

As in the BPS multi-centered attractors, solutions of Cf° are given by the 
single-centered non-BPS solutions — C°'°'{Q,z) with Q replaced by Qi. 
The solutions of twistors z and u are the same as the single-centered ones with 
charges Qnb replaces by Q^^. Among the aforementioned 2N—2 constraints, 
A'^ — 1 come from demanding that all centers have the same twistor z, which 
follows from the zero NUT charge condition at each center, and the other 



TV — 1 come from demanding that they have the same twistor u as in ( 89 1 . 
Solving these 2N — 2 constraints shows that all the charges {Qmbs} are the 
image of a single duality transformation /" on a multi-centered D4-D0 system 
{Qnbaoa}'- 

QNB,i — rQNB,40,i (92) 

The charges at different centers are all mutually local 

{QNB,^,QNB.J) ^0 (93) 

Like non-BPS single-centered attractors, the generic non-BPS multi-centered 
attractors cannot be generated via the "naive" harmonic function procedure, 
except for those with {QNB.i^h} being the image of a single F on a, pure 
D4-D0 system {Qnbao.i, hio}- 

tNB{x) ^ t*j,s{rQNBAO ^ E ^^^"^'^ + rfe40) (94) 

^ \X - Xi\ 



In summary, the non-BPS multi-centered attractors are drastically dif- 
ferent from their BPS counterparts: there is no constraint imposed on the 
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positions of the centers, but instead on the allowed charges QNB,i' they have 
to be mutually local. The result is that the centers can move freely, and there 
is no intrinsic angular momentum in the system. 

6 Conclusion and Discussion 

In this talk, we summarized the c;oustruction of generic single-centered and 
multi-centered extremal black hole solutions in theories whose 3D moduli 
spaces are symmetric coset spaces. In this construction, all attractors, both 
BPS and non-BPS, single-centered as well as multi-centered, are treated on an 
equal footing. The single-centered black hole attractors correspond to those 
null geodesies in Msd that are generated by exponentiating appropriate nilpo- 
tent elements in the coset algebra. The multi-centered black hole attractors are 
given by 3D solutions that live in certain null totally geodesic sub-manifolds 
of MzD- The construction of multi-centered attractors, even that of non-BPS 
ones, is merely a straightforward generalization of the single-centered con- 
struction. 

We presented a detailed computation in the theory of 4D A/^ = 2 supergrav- 
ity coupled to one vector- multiplet, whose 3D moduli space is the symmetric 
coset space G2(2)/'S'L(2,R)^. The attractor flow generators are third-degree 
nilpotent elements in the coset algebra. We expHcitly constructed generic at- 
tracrtor solutions, both single-centered and multi-centered, and showed that 
while the BPS attractors can be generated from the attractor moduli via the 
"naive" harmonic function procedure, the generic non-BPS attractors cannot 
be generated this way. 

In the n\/ = 1 model, besides the BPS generator, there is only one extra 
third-degree nilpotent orbit to serve as non-BPS flow generators. Hence there 
is only one type of non-BPS single-centered attractor. In models with bigger 
symmetric moduli spaces, there should be more than one type of non-BPS 
generator. These would give rise to different types of non-BPS attractor flows, 
which might have different stability properties. 

All multi-centered non-BPS attractors constructed in this work follow from 
the ansatz in which 3D gravity is assumed to decouple from the moduli. The 
multi-centered non-BPS black holes are found to be very different from their 
BPS counterparts: the charges of all centers are constrained to be mutually 
local, while the positions of centers are completely free. Thus the non-BPS 
multi-centered attractor is not a "bound state" and carries no intrinsic angular 
momentum. 

We would like to construct true multi-centered non-BPS "bound states" , 
i.e., solutions with constraints on the positions of centers but not on the 
charges. There are two possible ways to achieve this. First, one could adopt 
a more general ansatz in which 3D gravity is coupled to the moduli. For 

axisymmetric configurations, the inverse scattering method c;ould be used to 
perform an exact analysis. One could also search in models with bigger moduli 
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spaces. It is very likely that in bigger moduli spaces, there exist true multi- 
centered non-BPS "bound states" even within the ansatz with 3D gravity 
decoupled from moduli. We are also interested in the possibility of generating 
multi-centered non-BPS solutions with each center having different types of 
non-BPS generators k^B- 

Finally, with the hope of studying non-BPS extremal black holes in 4D 
M = 2 supergravity coupled to ny vector-multiplets with more generic pre- 
potcntial, we would like to generalize our method to non-symmetric homoge- 
neous spaces, and even to generic moduli spaces eventually. 
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